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We analyse numerically the critical behavior of a conserved lattice gas which was recently introduced 
as an example of the new universality class of absorbing phase transitions with a conserved field 
[Phys. Rev. Lett. 85, 1803 (2000)]. We determine the critical exponent of the order parameter as 
well as the critical exponent of the order parameter fluctuations in _D = 2, 3, 4, 5 dimensions. A 
comparison of our results and those obtained from a mean-field approach and a field theory suggests 
that the upper critical dimension of the absorbing phase transition is four. 



I. INTRODUCTION 

The scaling behavior of directed percolation is recog- 
nized as the paradigmatic example of the critical behav- 
ior of several non-equilibrium systems which exhibits a 
continuous phase transition from an active state to an 
absorbing non-active state (see for instance [|l],D). Such 
systems are common in physics, biology, as well as cat- 
alytic chemical reactions. This widespread occurrence 
corresponds to the well known universality hypothesis 
of Janssen and Grassberger that models which exhibit a 
continuous phase transition to a single absorbing state 
generally belong to the universality class of directed per- 
colation 

Recently Rossi et at. introduced a conserved lattice gas 
(CLG) with a stochastic short range interaction that ex- 
hibits a continuous phase transition to an absorbing state 
at a critical value of the particle density The CLG 
model is expected to belong to a new universality class of 
absorbing phase transitions characterized by a conserved 
field. Similar to the above hypothesis the authors conjec- 
tured that " all stochastic models with an infinite number 
of absorbing states in which the order parameter evolu- 
tion is coupled to a nondiffusive conserved field define a 
unique universality class" ||]. Besides the CLG model 
the authors considered the conserved threshold transfer 
process model as well as a modification of the stochastic 
sandpile model of Manna Q and observed numerically 
compatible values of the critical exponents. Furthermore 
a reaction-diffusion model was introduced in Q which is 
expected to belong to the same universality class. This 
reaction diffusion model allows to derive a field theoreti- 
cal description which is expected to represent the critical 
behavior of the whole universality class |^]. 

In this work we consider for the first time the scal- 
ing behavior of the CLG model and therefore the critical 
behavior of the new universality class in higher dimen- 
sions. We determine numerically the critical exponent of 
the order parameter as well as the exponent of the order 
parameter fluctuations. Our results show that the val- 
ues of the exponents depend on the dimension for D < A. 



Above this dimension we observe a mean-field scaling be- 
havior. Thus our results suggest that the upper critical 
dimension of the CLG model is four as already predicted 
from the field theoretical approach ||^. 

II. D=2 

We consider the CLG model on Z?-dimensional cubic 
lattices of linear size L. Initially one distributes randomly 
iV = pL particles on the system where p denotes the par- 
ticle density. In order to mimic a repulsive interaction a 
given particle is considered as active if at least one of its 
2D neighboring sites on the cubic lattice is occupied by 
another particle. If all neighboring sites are empty the 
particle remains inactive. Active particles are moved in 
the next update step to one of their empty nearest neigh- 
bor sites, selected at random. Starting from a random 
distribution of particles the system reach after a tran- 
sient regime a steady state which is characterized by the 
density of active sites p^. The density pa is the order pa- 
rameter of the absorbing phase transition, i.e., it vanishes 
if the control parameter p is lower than the critical value 
Pc. In contrast to the work of Rossi et al. who used a 
parallel update scheme, we applied in our simulations a 
random sequential update, i.e., all active sites are listed, 
and then updated in a randomly chosen sequence. 

We consider in two dimensions simple cubic systems of 
linear size L — 32, 64, 128, . . . , 2048 with periodic bound- 
ary conditions. Starting with a random configuration 
of particles, a sufficient number of update steps has to 
be performed to reach the steady state where the num- 
ber of active sites fluctuates around an average value (see 
Fig. Approaching the transition point, more and more 
update steps are needed to reach this steady state. For 
instance in the case L = 2048 we use 2 10® update steps to 
"equilibrate" the system. In the steady state the number 
of active sites is monitored for 5 10*^ update steps. This 
procedure is repeated in all dimensions for at least 10 
different initial configurations. From this data we deter- 
mine the average density of active sites (pa) as well as its 
fluctuations 
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FIG. 1. The density of active sites pa as a function of time 
(number of update steps) for a certain value of p. After a 
transient regime, which depends on the initial configuration, 
the density of active sites fluctuates around the steady state 
value (pa) (dashed line). 



(1) 



In Fig. ^ we plot the average density of active sites 
as a function of the particle density p for various system 
sizes L. As one can see (pa) tends to zero in the vicinity 
of p « 0.345. Assuming that the scaling behavior of the 
density of active sites is given by 



(Pa) - {p-PcY 



(2) 



one varies pc until one gets a straight line in a log-log 
plot. Convincing results are obtained for p^ = 0.34494 ± 
0.00003 and the corresponding curve is shown in Fig.|^. 
For Pc = 0.34491 and pc — 0.34497 we observe significant 
curvatures in the log-log plot (see inset of Fig. ||). In 
this way we estimate the error-bars in the determination 
of the critical density. A regression analysis yields the 
value of the order parameter exponent (3 — 0.637 ±0.009. 
Rossi et al. reported the values (3 — 0.63 ± 0.01 and 
Pc = 0.28875, obtained from simulations with a parallel 
update scheme and of smaller system sizes [L < 512) 
Thus we see that the different update scheme affects only 
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FIG. 2. The average density of active sites (pa) as a func- 
tion of the global particle density p for various system sizes L 
m D = 2 dimensions. 
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FIG. 3. The average density of active sites (pa) as a 
function of p — pc in D = 2. The symbols mark dif- 
ferent system sizes L (see Fig. ^) . The dashed line corre- 
sponds to a power-law fit with pc = 0.34494 ± 0.00003 and 
13 = 0.637 ± 0.009. In the inset we display the same data 
for Pc = 0.34491 and pc = 0.34497. Compared to the above 
value (solid line) both curves are characterized by significant 
curvatures in the plotted log-log diagram. For the sake of 
simplicity we plot in the inset lines instead of symbols. 



the value of the critical density but not the critical ex- 
ponent (3. Furthermore the value of (3 differs from the 
corresponding value of the directed percolation univer- 
sality class /3dp = 0.584 ± 0.004 (see [§). 

The fluctuations of the order parameter [Eq. (|l])] are 
plotted in Fig.^. Approaching the transition point Apa 
increases and diverges at pc. Close to the critical point 
the fluctuations scale as 



Apa 



(3) 



in the active phase (see inset of Fig. ^). Using a regression 
analysis one gets 7 = 0.384 ± 0.023. 



3.0 t- < 

: i 
r 

2.0 t 



0=2 



10 10 



0.34 0.35 0.36 0.37 0.38 0.39 

P 

FIG. 4. The fluctuations of the order parameter Apa as 
a function of the global density p in D — 2. The symbols 
mark different system sizes (see Fig. ^) and the dashed line 
corresponds to the critical density pc = 0.34494. The inset 
displays the fluctuations as a function of p — p^ . The dashed 
line corresponds to a power-law behavior [Eq. (^)] with an 
exponent 7 = 0.384 ± 0.023. 
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FIG. 5. The average density of active sites (pa) as a func- 
tion of p — pc for various system sizes L in D = 3. The dashed 
hne corresponds to a power-law fit with pc = 0.2179 ± 0.0001 
and /3 = 0.837 ± 0.015. 



universality class /3dp = 0.81 ± 0.01 (see 

Similar to the two dimensional case the order parame- 
ter fluctuations Apa display a maximum at the transition 
point but the dependence of Apa on the density of par- 
ticles is not clear. In Fig.^ we plot the fluctuations in 
a log-log plot. It seems that the asymptotic behavior 
corresponds to a power-law [Eq. (|^)] with an exponent 
7 — 0.18 ±0.06. But as the inset of Fig. ^ shows the data 
are also consistent with the assumption that the fluc- 
tuations are characterized by a logarithmic divergence 
(i.e. 7 = 0) at pc according to 



Apa 



In (p - Pc 



(4) 



Thus the numerical data indicate that the fluctuations 
Apa diverge at the critical point but the data can be 
interpreted either as a power-law with a small exponent 
or as a logarithmic growth. Further investigations are 
needed to clarify this point. 



III. D=3 



IV. D=4 



For the three dimensional model system sizes from 
L = 16 to L = 160 are considered. Close to the transi- 
tion point 5 10^ update steps are used to reach the steady 
state in the case of the largest system size. The ob- 
tained results for the density of active sites (pa) are shown 
in Fig. ||. A straight line in a log-log plot is obtained 
for Pc ~ 0.2179 ± 0.0001 and the corresponding value 
of the order parameter exponent is /? = 0.837 ± 0.015. 
The error-bars of the critical density are determined in 
the same way as for = 2. The value of the order 
parameter exponent agrees within the error-bars with 
the corresponding values of the Manna sandpile model 
{/3 = 0.84 ± 0.02) and of a reaction diffusion model 
(/3 = 0.86 ± 0.02) which are expected to belong to the 
same universality class Q. Our result differs slightly 
from the corresponding value of the directed percolation 



In order to analyze the scaling behavior of the four 
dimensional CLG model we performed numerical simula- 
tions with system sizes L G {8, 16, 32, 48}. In the case of 
the largest system size 6 10^ updates steps were used to 
reach the steady state. Plotting the values of the aver- 
age particle density (pa) in a log-log plot no straight line, 
i.e., no pure power-law behavior could be observed. To 
illustrate this behavior we plot in Fig. ^ the logarithmic 
derivative 



d In (pa) 
9 In (p - Pc) ' 



(5) 



which can be interpreted as an effective exponent PcS- If 
the scaling behavior of the active site density is given by 
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FIG. 6. The fluctuations of the order parameter Apa as a 
function of p — pc in D = 3. The symbols mark different 
system sizes (see Fig.^. The data can be interpreted either 
as a power-law with a small exponent 7 (the dashed line cor- 
responds to 7 = 0.18 ± 0.06) or as a logarithmic growth [see 
Eq. (^) and inset, where the dashed line is to guide the eye]. 
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FIG. 7. The logarithmic derivative [Eq. ^] as a function 
of p — pc. The logarithmic derivative can be interpreted 
as an effective exponent PeS- The figure shows that the 
four-dimensional exponent does not display a saturation as 
the exponents of the three and five dimensional model do. 
The dashed lines correspond to the three dimensional value 
P = 0.837 ± 0.015 and the mean-field value /3 = 1, respec- 
tively. The shadowed region marks the uncertainty of the 
determination of j3. 
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FIG. 8. The density of active sites in D — i 
rescaled according to Eq. ^ . The assumed asymptotic 
behavior (straight line in a log-log plot) is obtained for 
= 0.1571 ± 0.0002 and B = 0.39. Thus the scaling be- 
havior of the order parameter of the four dimensional model 
is governed by the mean-field exponent (3—1 modified by 
logarithmic corrections. The dashed line is just to guide the 
eye. 
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FIG. 9. The fiuctuations of the order parameter Apa as 
a function of {p — pc)| In (p — pc)\^ in D = 4. [see Eq. (^)]. 
Nearly straight lines are obtained for F = 1.66 ± 0.1. Thus 
the scaling behavior of the order parameter of the four dimen- 
sional model is governed by the mean-field exponent 7 = 
modified by logarithmic corrections. The dashed line corre- 
sponds to a linear fit with the slope —0.643 and the intercept 



Ap!^ 



0.325. 



Eq. (g) the logarithmic derivative tends to the value of /3 
for p — Pc ^0. This behavior is observed in the three 
dimensional case (see Fig. . For D = A the logarithmic 
derivative displays no saturation for p —^ pc, i.e., the 
scaling behavior of the four dimensional model can not 
be described by a simple power-law behavior. Significant 
corrections to the usual scaling behavior [Eq. (y)] occur 
for instance at the upper critical dimensional where the 
scaling behavior is governed by the mean-field exponents 
modified by logarithmic corrections. 

Recently a modified version of the CLG model was 
introduced where the active particles are distributed to 
randomly chosen empty lattice sites Since the ran- 
domness of the particle hopping breaks long range spatial 
correlations this model is expected to be characterized by 
the mean-field scaling behavior of the CLG model. Map- 
ping the dynamics of this random hopping CLG model 
to a simple branching process one can derive the critical 
exponent (3=1. This value of the order parameter expo- 
nent is also obtained from a field theoretical description 
of the CLG model ^ and was already predicted from 
a phenomenological field theory [|l^ of the so-called fix- 
energy stochastic sandpile model which is expected to be 
in the same universality class. 

Thus we assume that the scaling behavior of the order 
parameter is given in leading order by the ansatz 

(Pa) ^ (P-Pc)''|ln(p-Pc)r (6) 

with (3 = 1. Therefore we varied in our analysis B and 
Pc until we get the expected asymptotic behavior. The 
best result is obtained for B = 0.39 and pc = 0.1571 ± 
0.0002 and the corresponding scaling plot is shown in 
Fig. ||. As one can see our data are consistent with the 
assumption that the asymptotic scaling behavior of the 
four dimensional model obeys Eq. (^). 

The mean-field behavior of the fluctuations is charac- 
terized by 7 = which corresponds to a jump [^ . Taking 



the logarithmic corrections into account we assume that 
the asymptotic scaling behavior of the fluctuations obeys 
the ansatz 

Apa = Api"' - const [p - Pc) I In (p - p^)!^. (7) 

Plotting the fluctuations as a function of (p — 
Pc) I In (p — pc)|^ one varies the correction exponent F 
until one gets a straight line. Convincing results are ob- 
served for F — 1.66 ± 0.1 and the corresponding plot 
is shown in Fig.^ (for F > 1.76 and F < 1.56 we ob- 
serve significant curvatures). The extrapolation to the 
vertical axis yield the fluctuations at the transition point 
Ap™ = 0.325 ± 0.005. 

Thus we get that the scaling behavior of the four di- 
mensional CLG model is characterized by the mean-field 
exponents modified by logarithmic corrections and we 
conclude therefore that the upper critical dimension of 
the CLG model is Dc — 4. This value agrees with the 
conjecture of a field theory 

V. D=5 

In the case of the five dimensional model we considered 
system sizes from L = 4 up to L = 18. In the latter case 
10^ update steps were used to reach the steady state close 
to the transition point. The obtained values of the order 
parameter are plotted in Fig. The average density 
of active sites seems to vanish linearly at the transition 
point. This is supported by the logarithmic derivative 
[Eq. (^)] which is plotted in Fig. ^ The effective exponent 
saturates in the vicinity of the mean-field solution (3 = 
1. One has to admit that the small available system 
sizes in d = 5 (L < 18) hinder an analysis of the pure 
critical scaling behavior. But as we will see the data are 
sufficient to indicate that the asymptotic scaling behavior 



Physical Review E 64, 016123 (2001) 



5 



■ L=4 
1 L= 6 

■ L= 8 

■ L=10 
1 L=12 
' L=14 

L=16 
L=18 



0.00 I 

0.10 



D=5 



0- 10 1 

0.25 0.30 



TABLE I. The critical density pc and the critical exponents 
(3 and 7 of the CLG model for various dimensions D. The 
symbol * denotes logarithmic corrections to the power-law 
behavior. In the case of the three dimensional model the data 
of the fluctuation could be interpreted as a small exponent or 
as a logarithmic growth (see text). 
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3 0.2179 ± 0.0001 
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FIG. 10. The average density of active sites (pa) as a func- 
tion of the density p for various system sizes L in D = 5. The 
inset displays (pa) /(p — pc) as a function of p — pc. The satu- 
ration for p — Pc indicates that the asymptotic behavior 
of the order parameter agrees with the mean field behavior 
(pa) ~ p - Pc. The critical density is pc = 0.1230 ± 0.0004. 



is given by /3 = 1. Therefore, we plot in the inset of 
Fig. |l0| (pa)/(/0 ~ Pc) as a function of p — pc- Approaching 
the transition point the rescaled density of active sites 
indeed saturates. Thus the five dimensional CLG model 
is characterized by the mean-field scaling behavior 



(Pa) 



P- Pc- 



(8) 



The fluctuations of the order parameter Apa are plot- 
ted in Fig.|ll|. As one can see the fluctuations are char- 
acterized by a jump at the transition point. Following 
the mean-field behavior one expects that the asymptotic 
behavior of the fluctuations is given by 

Apa = Ap!,"' - const (p-pe). (9) 

In order to confirm this ansatz we plot Apa/ (p — Pc) as a 



function of p — p^. If the above ansatz is valid the corre- 
sponding curves have to display an asymptotic power-law 
behavior with the exponent —1 and the prefactor Ap^^"'. 
This behavior is indeed observed (see inset of Fig.|ll|) and 
we get Ap(f' = 0.231 ± 0.006. 

Thus the asymptotic behavior of the five dimensional 
CLG model is characterized by the mean-field exponents. 
This behavior strongly supports the above conclusion 
that the upper critical dimension of the CLG model is 
four. 



VI. CONCLUSIONS 

We analyse numerically the critical behavior of a CLG 
model in various dimensions. The values of the criti- 
cal density and of the critical exponents are determined 
and the results are listed in Table |. Our analysis sug- 
gests that four is the upper critical dimension of the CLG 
model, i.e., the critical exponents depend on the value of 
the dimension for _D < 4. Above this value we observe a 
mean-field like scaling behavior. 
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FIG. 11. The fluctuations of the order parameter Apa as a 
function of the density p in D = 5. The symbols mark differ- 
ent system sizes L (see Fig.|l^. The dashed line marks the 
value of the critical density. The inset displays Apa/(p — pc) 
as a function of p — pc. The asymptotic behavior is character- 
ized by a simple (p — Pc)~^ dependence (dashed line) which 
corresponds to the mean- field behavior (see text). 
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